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ABSTRACT 

The observed luminosity function (LF) of satellite galaxies shows several interesting features that require a 
better understanding of gas-thermodynamic processes and feedback effects related to reionization and galaxy 
formation. In galaxy clusters, the abundance of dwarf galaxies is in good agreement with the expectation based 
on the subhalo mass function, whereas in galaxy groups, the relatively small abundance of dwarfs conflicts with 
theoretical expectations. In all halo systems, there is a dip in the abundance of galaxies with luminosities in the 
range ~2x 10 8 L Q to 10 10 L Q , corresponding to subhalo mass scales between ~5x 10 10 M Q to few times 
10 11 Mq. Photoionization from reionization has been used to explain statistics of the dwarf population, with 
larger systems forming prior to, and smaller systems forming subsequent to, reionization. The observed dip in 
the LF is an imprint of small dwarf galaxies (<~ 2 x 10 8 L Q ) that formed prior to reionization. The galactic 
winds powered by supernovae in these dwarf galaxies propagate energy and metals to large distances such that 
the intergalactic medium is uniformly enriched to a level of 10~ 3 Z Q , as observed in the low-redshift Lya forest. 
The associated energy related to this metallicity raises the intergalactic medium temperature and the Jeans mass 
to a range 10 10 - 10 11 M Q at z ~ 3.4-6.0. Because the epoch of nonlinearity for halos in this mass range is at 
Z > 3 A- 4.4, their gas content, hence star formation, is greatly suppressed on average and leads to a dip in the 
observed LF at z = 0. Larger halos (M > 10 11 M Q ), becoming nonlinear at z < 3.4-4.4, have masses topping 
the Jeans mass, where subhalo mass function based LF is again in agreement with observations. 
Subject headings: cosmology: observations — cosmology: theory — galaxies: clusters: general — galaxies: 
formation — galaxies: fundamental parameters 



1. INTRODUCTION 

The complex physical processes associated with reioniza- 
tion of the intergalactic medium (IGM) at a redshift above 6 is 
expected to leave characteristic scales and features in the star 
formation h istory (e.g..lCenll2003l) . in the supernovae distri- 
bution (e.g., Mesinger et al. 2005), and, potentially, in the lu- 
minosity distributions of galaxies. In the case of galaxy statis- 
tics, for example, the feedback related to supernovae heating 
in small mass halos and photoionization during reionization 
llBenson et all2002 aS). have been used to explain the flattened 
faint-end slope of the galaxy luminosity function (LF). 

Since one averages galaxy statistics over large volumes and, 
thereby, averages over any inhomogeneities and differences 
in time scales and dispersions coming from galaxy formation 
and evolution processes, it is not possible to address detailed 
physics related to reionization with the galaxy LF alone. On 
the other hand, the luminosity distribution of satellites in dark 
matter halos, as a function of the halo mass, may be an ideal 
probe of reionization physics. In this respect, the lack of an 
abundant population of low luminosity galaxies in the local 
group, relative to expectations from cold dark matter cos- 
mological mod els, has been explained in terms of photoion- 
ization JBullock et alJl2000HBenson et al]l2003l The subse- 
quent squelching of galaxy formation in dark matter halos be- 
low a certain mass scale, corresponding to the temperature to 
which IGM is heated, explains the environmental dependence 
of the faint-end slope of the cluster LF (Tullv et al. 2002). 

In addition to effects related to photoionization, the 
cluster satellite LF should also show signatures of feed- 
back associated with starformation. The satellite luminos- 
ity distribution can be described thro ugh the halo occu- 
pation number dCoorav & Shethl 120021) . minus the central 



galaxy, conditioned in terms of satellite luminosity. Here, 
we construct an empirical mod el for the conditional lu- 
minosity function (CLF; Yang, Mo, & van den BoscU 120031 
Coorav & Milosavlievic 2005: lCooravI f2005) of satellites in 
dark matter halos, based on the subhalo mass function and 
compare to observed measure ments; at th e brig ht-end, we 
make use of CLFs m easured by Yang et al. (2005) using the 
2dFGRS llCole et alJl2001l) galaxy group catalog, and extend 
thi s comparison to the faint-end usi ng cluster LFs me asured 
by Trentham & Hodgkh] J2001 and lHilker et all J200I . We 
argue that cluster LFs show two scales, one associated with 
photoionization at halo mass scales around 5 x 10 10 M Q , re- 
sulting in the disappearance of dwarfs satellites in galaxy 
groups relative to clusters, and another scale related to an 
overall suppression of galaxy formation in subhalos below 
10 11 Mq, independent of the total system mass. 

This Letter is organized as follows: In § |2j we describe 
the construction of the satellite CLF of dark matter halos. In 
§|3j we compare our LF with the observed LFs of satellites in 
groups and clusters and discuss interesting physics that could 
explain the observed features. We make use of the current 
concordance cosmological model ( Sper geTet alJ20 03). 

2. SATELLITE LUMINOSITY FUNCTION 

The CLF, denoted by <1>(L|M), is the average number of 
galaxies with luminositie s between L and L+dL that reside in 
halos of mass M ( Yang. Mo. & van den Bosch oorav 
120051) . Following Coorav & Milosavlievic (2005), we sepa- 
rate the CLF into terms associated with central and satellite 
galaxies, <1>(L| M) = $ C (L|M) + $ S (L|M). In previous studies 
(Coorav 2005), central galaxy CLF was described with a log- 
normal distribution in luminosity with a mean L C (M) and dis- 
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FIG. 1 . — The CLF of dark matter halos as a function of the halo mass. From (a) to (d) we show four mass ranges from group to cluster scales, with mass ranges 
(in logarithmic values) labeled within square brackets in each of the figures. The circles at the bright-end correspond to the measurements in the foj-band by 
Yang et al. 1 2005) using the 2dFGRS i Cole et al. 2001) group catalog with a division to satellites (open circles) and satellites and central galaxies (solid circles). 
The long-dashed line is the central galaxy CLF following Coorav 1 2005). The top solid line shows the expected CLF of satellites based on the subhalos mass 
function. The lower solid line is the satellite CLF when the subhalo mass function has a turnover at masses below 10" Mq. The faint-end data comes from 
measurements in the B-band by Trentham & Hodskin 1 2002) for Ursa Major, in (a) with a total mass around ~ 8 X 10 12 Mq, Virgo, in (b) with a total mass of 
~ lO 14 Mq, and Coma, in (c), where we take the total mass to be between (7 and 10) X 10 14 Mq. In (a) and (b), we also show the faint-end LF of Fornax (from 
I). In published literature, the total mass of Fornax varies between (4.3 to 8.1) X 10 12 Mq, based on X-ray data i Jones et al. 1997), to ~ 7 X 10 1 



|Hilke^^y2003). In published literature, the total mass of Fornax varies between (4.3 to 8.1) X 10 12 Mq, based on X-ray data i Jones et al. 1997), to ~ 7 X 10 13 
(e.g., Bekki et al. 2003). If the low halo mass is correct for the Fornax group, then Fornax formed earlier than Ursa Major and before reionization was complete. 
If the high mass is correct, given the presence of some dwarf galaxies, but below the level expected, Fornax probably was forming during reionization. The data 
can also be used to argue that the reionization process had inhomogeneities at size scales below that collapsed to form Fornax. Note the difference in the y-axis 
scale between the two top panels and the two bottom panels. 



persion S c , while the satellite CLF is assumed to be a power 
law. Here, we focus on the satellite CLF and, instead of an a 
priori assumption on a power-law CLF, we model it using the 
subhalo mass function. 

In this approach, each satellite in a subhalo mass M s has a 
log-normal luminosity distribution of 
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where the normalization <j)(M s ) is fixed such that 
J <&{L\M s )LdL = L S (M S ). Given the luminosity distribu- 
tion of each satellite, the CLF of satellites, in a parent halo 



mass M, is 



$ S (L|M) : 



dn s (M s \M) 
(L\M S ) ' dM s , 



o dM. 



(2) 



where dn s (M s \M) jdM s is the subhalo mass function of dark 
matter halos given the parent halo mass M. Here, we use the 
analytical form 

dn s {M s \M) _ 7 f M A~ a (Jh}\ m 

dM s P 2 MT(2-a)\l3M) 6XP I /3M J ' ( > 



where a = 1.91, /3 = 0.39, and 7 = 0.18 dVale & Ostrikerl 
2004). Our concl usions do not change signifi cantly if we use 
the description of Ivan den Bosch et alJ ( 120051) where a and 7 
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are functions of mass, M, with a varying roughly over 10% 
as M is varied from group to cluster mass scales. Since CLF 
measurements are average d over a sample o f dark matter halos 
in a narrow range in mass dYang et al.l2005l) . we also calculate 
the mass-averaged satellite CLF by averaging over the dark 
matter halo mass function, dn/dM (Sheth & Tormen 1999), 
over the mass range of interest. 

In our model, an important ingredient is the L S (M S ) rela- 
tion which describes the luminosity of a subhalo given the 
subhalo mass. Here, we fo llow the procedure related to 
modeling the field galaxy LF (Coorav & Milosavlievic 2005), 
and employ a fitting function to describe the relation be- 
tween the lumin osity of a halo and the dark matter mass 
of that halo. In Vale & Ostriker (2004), this relation was 
obtained through a model description of the 2dFGRS LF 
l lNorberg et all20 01 ) using the global subhalo mass function, 
n S h(M s ) = f dn/dM dn s (M s \M)/dM s dM. The relation is 



L(M) = Lq 



(M/Mif 



[b+{M/M x y d yi d 



(4) 



The relevant parameters for the fo/-band, as appropriate for 
2dFGRS data, are L = 5.7 x 10 9 L©, M] = 10"Afo, a ; = 4.0 , 
b = 0.57, c = 3.72, and d = 0. 23 (IVale & Ostrikerl l200l . 
Though this relation was used in Coorav (2005) to describe 
the field galaxy LF, whose statistics are dominated by central 
galaxies, the same relation should also remain valid for sub- 
halos as well. The remaining parameter in our model for the 
satellite LF is £ and we set this to be 0.17 based on the value 
nee ded to explain the exponential drop-off in the field galaxy 
LF (ICoora\l2005l) . 

Note that at low subhalo masses, dn s /dM s oc M~ a where 
a = 1.91, independent of the halo mass, or varies from 2.0 
to 1.9 wh en parent halo mass varies from M ~ 1O U M0 to 
10 15 M Q (van d en Bosch et al.1 120051) . The satellite CLF is 
$ S (L|M) ~ / 4> s (L\M s )(dn s /dM s )dM s . If L s oc M%, we can 
write $ S (L|M) - f L'~ l ~ a / T > +l / T > 5(L' - L)dL' , where we have 
ignored the scatter in the L s -M s relation by setting £ — > 0. 
The faint-end of the satellite LF then scales as $ S (L|M) oc 

T—l—a j r]+\ / rj 

3. RESULTS AND DISCUSSION 

In Figure 1, we show b r-band CLFs based on the 2dFGRS 
group catalog (Yang et al. 2005) and an extension to the faint- 
end based on nearby cluster LFs in the B-band. The satel- 
lite CLFs show several interesting trends. As one moves to a 
higher mass for the central halo, from (a) to (d) in Figure 1, 
one finds the faint-end of the satellite CLF to be filled with 
dwarf galaxies. Since the 2dFGRS is incomplete below ab- 
solute magnitude s of Mh, of -17, we make use of B-band LF 
measurements by Trentham & Hodgkin (2002), whose faint- 
end statistics are generally dominated by the dwarf galaxy 
population. At the high mass end of halos, corresponding to 
clusters like Virgo, with an assumed halo mass of 10 14 M Q , 
and Coma, with mass between ~ (7 and 10) x 10 14 M©, satel- 
lites begin to trace the expected slope predicted by the subhalo 
mass function. On the other hand, in low mass groups, such as 
Ursa Major, with an assumed mass of (5 to 8) x 10 12 M©, the 
faint-end dwarf population does not trace the subhalo mass 
function. The statistics of t he faint-end dwarf galaxy popula- 
tion have been discussed in Tullv et al. (2002) in the context 
of photoionization effect resulting from reionization. These 
authors argue that massive clusters such as Virgo started to 
form prior to reionization, while low mass groups such as 
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FIG. 2. — The subhalo mass function assuming a parent halo mass of 10 
Mq . The solid line show the expectation based on CDM calculations (equa- 
tion 0. The middle line shows the flattening needed to suppress the pres- 
ence of galaxies fainter than 3 X 10 9 Lq or the corresponding mass scale 
of roughly 10" Mq. The flattening is achieved through a smooth function 
discussed in Section 3. The third line labeled "sub-halos with dwarfs" is the 
subhalo mass function necessary to explain the luminosity distribution at the 
faint-end of L < 3 X 10 8 Lq when the parent halo mass is above 10 14 Mq. 



Ursa Major, where the abundance of dwarf galaxies is rela- 
tively smaller, formed subsequent to reionization. As shown 
in Figure 1(c) and (d), the faint-end dwarf LF is associated 
with subhalo masses below 5 x 10 10 M©. The photoionized, 
and heated, gas is not expected to cool in halos of below this 
mass scale after reionization. Such an argument is consistent 
with Figure 1 . 

Based on satellite LFs alone, we find that reionization may 
be more complex than simply assuming that the universe 
reionized completely at a single redshift. For example, in 
Figure 1(a) and (b), we also show the LF of dwarf galaxies 
in Fornax. We plot these data in both panels due to a large 
uncertainty, or variation, in the quoted total mass of Fornax 
in the literature. Regardless of the exact mass of Fornax, 
the presence of more dwarfs than Ursa Major suggests that 
Fornax formed prior to the latter system and probably dur- 
ing reionization such that gas managed to cool in a fraction 
of small dark matter halos while the rest was affected. The 
large scatter in dwarf population of similar mass groups may 
be evidence for inhomogeneous reionization or non-uniform 
feedback processes. 

What is new and intriguing is this: when compared to 
the expectation based on th e sub-halo mass function and the 
luminosity-mass relation of Vale & Ostriker (2004), one sees 
a relative decrease in the number of satellite galaxies at lumi- 
nosities below ~ 3 x 10 9 L©, corresponding to mass scales 
below few times 10 11 M©, in all halo systems. Since at these 
mass scales, 3 < rj < 4, we expect $ S (L|M) to scale as L" L2 
to L" 1 3 ; we do not expect the faint-end slope to be flatter than 
L" 1 2 , unless the subhalo mass function slope is changed; for 
$ S (L|M) to be luminosity independent, a — ► 1. Thus, to sup- 
press galaxies at luminosities below 3 x 10 9 L©, we include 
an efficiency function to the subhalo mass function to char- 
acterize the subhalo mass distribution where satellite galaxies 
present. Here, we take an analytical description of the form 
f(M s ) = 0.5(1 +erf[(logM s -logM f .)]/cr), such that f(M s ) -> 
when M s > M c and f(M s ) — > when M s < M c . To explain 
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the flattening of the satellite CLF, we set M c = 10 11 M and 
a = 1.5; with such a broad dispersion, f(M s ) dn{M s \M)/dM s 
flattens (see, Figure 2) instead of becoming zero even when 
M ~ 10 9 M Q ; the flattening of the satellite CLF does not im- 
ply that all subhalos below some critical mass scale is af- 
fected, but one sees a broad distribution of subhalo masses 
where, statistically, satellite galaxies are not present. Instead 
of modifying the subhalo mass function, we can also vary the 
L S (M S ) relation and set a steep slope for r\. While in this case 
all subhalos contain galaxies, the average luminosity would 
be lower for subhalos in groups relative to same mass subha- 
los in clusters. In such a scenario, it is also hard to understand 
the sudden appearance of dwarfs in clusters, given the dip in 
the LF of giants regardless of system mass. 

We now offer a physical explanation for this unique feature. 
Supernovae-powered winds from abundant dwarf galaxies at 
L < 2 x 10 s M are expected to be strong (Dekel & Silk 1986; 
Mori et al. 2002). We adopt the view that these winds trans- 
port metals and energy into the IGM. Following Cen & Bryan 
(2001), the temperature of the IGM is 



T SN (z)= 1.3 x 10 4 K 
M c 



0.2 M, 



0.3 



1.2 x 10 51 erg 
[C/H] \ ( 4 



1 x 10" 



1+z 



(5) 



where Mc is the mass of carbon ejected by one supernova; E$n 
is the total energy output of one supernova; r\ is the fraction of 
that energy that is eventually deposited in the IGM in the form 
of thermal energy (Mori et al. 2002), and [C/H] is the ratio 
of carbon number density to hydrogen number density of the 
gas in solar units. We have assumed that energy is deposited 
at some high redshift, perhaps z ~ 6 - 9, and the temperature 
of the IGM subsequently decays adiabatically. We adopt all 
the fiducial values for Esn, Mc, t] and [C/H], which gives us 
an added temperature to the IGM of 1.3 x 10 4 K at z ~ 3; this 
explains the Doppler width issue of the Lya forest (Cen & 
Bryan 2001). 

Figure 3 shows the evolution of the nonlinear mass in the 
standard cold dark matter model (solid curve), the Jeans mass 
for two extreme conditions (two short dashed curves), and the 
temperature of the IGM (long dashed curve). We see that, 
depending on detailed gasdynamics, the majority of halos 
formed at z> 3.4-4.4 and with masses 10 1() -6x 10 I0 /T 1 M Q , 
are significantly deprived of gas, potentially explaining the 
dip in the satellite LFs seen in Figure 1 . 

Our results suggest that most subhalos of 10 10 M Q to 10 11 
Mq in massive clusters should not have galaxies in them 
(so-called "dark halos"). Lensing-based studies indicate the 



agreement between dark matter subhalo mass function and the 
mass function of subhalo s associated with galaxies in clusters 
whenM. s > 3 x 10 11 M Q JNataraian & Springell2004l) . While 
the observed difference below this mass scale is considered 
to be due to an observational limitation, it is useful to extend 
lensing studies to measure the subhalo mass function below 
10 11 Mq. The dip in the LF could be detected as a dip in the 
mass function of subhalos with galaxies (Figure 2). Further- 
more, a large sample of satellite LFs, down to dwarf galaxy 
magnitudes and with better determined parent halo masses, 
have the potential to address more detailed physics of reion- 
ization and galaxy formation processes through a better char- 
acterization of features in the CLF. 



AC thanks Frank van den Bosch for helpful correspondence 
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FIG. 3. — The solid line shows the nonlinear mass scale in the standard 
LCDM model as a function of redshift. The two short-dashed curves show 
Jeans mass assuming dark matter is uniformly distributed (upper curve) and 
for gas at the virial density (i.e., 200 times the mean density) assuming gas 
and dark matter cluster in the same way (lower curve). The long-dashed curve 
shows the evolution of IGM temperature (values labeled on the right y-axis). 
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